NIKHEF-98-021 

Some problems in loop calculations 



J. A.M. Vermaseren 

NIKHEF, P.O. Box 41882, 
1009 DB, Amsterdam, The Netherlands 

We discuss some of the problems that may occur in the calculation of complicated Feynman 
diagrams. These include the group independent evaluation of color factors, and the summation 
techniques that are needed for the expansion of diagrams into their Mellin moments. 



1. Introduction 

The ever increasing accuracy of the high energy experiments forces theorists to calculate 
perturbative quantities to higher and higher order. For high order graphs there are however 
not only the integrals to worry about. One of the problems is the organization of calculations 
with astronomical numbers of diagrams. This we will not discuss here. Another problem 
concerns the color factors. Traditionally they have been evaluated for QCD with an algorithm 
that was specific for the SU(N) groups Q Q. These results could then be rewritten in terms 
of Casimir invariants and applied to other groups as well. This procedure works well as long 
as the only relevant invariants are Cf and Ca- Recent calculations (e.g. ref. [||) however 
went beyond this and needed additional invariants. Hence new algorithms were needed. We 
will discuss them here. 

Another problem that arises concerns the hadron structure functions in deep inelastic 
scattering. Here the structure functions can be computed to two loops but for a complete 
NNLO analysis one will also need the three loop anomalous dimensions. Of these anomalous 
dimensions only some Mellin moments are known thus far. This raises several points simulta- 
neously. The first is the computation of the Mellin moments, once the exact result in x-space 
is known. The second problem is that one would like to use these moments for at least a 
partial analysis. And thirdly one may wonder whether it could be possible to evaluate all mo- 
ments simultaneously. Or in other words: can one compute these diagrams in Mellin space? 
This technique has been used in the past for several calculations H || 0, but each time the 
calculations were done mainly by hand, in which case one can use many tricks to come to an 
answer. Additionally these techniques have never been used to the level of complexity one 
needs for the anomalous dimensions in NNLO. In the case of large numbers of diagrams one 
will have to solve the problem rather thoroughly. We will discuss some of this and show some 
progress in at least the first point. 

2. Color factors 

The evaluation of color factors is a purely mathematical problem. Yet most of the lit- 
erature about it has been written by physicists. The most widely used paper is the one by 



2 



Cvitanovic BJ. In this paper explicit properties of the fundamental and adjoint representa- 
tions of various groups are used to break the traces down and obtain expressions for the traces 
in terms of N (for SU(N), SO(N) and Sp(N)) or just numbers (for G 2 , F 4 , E 6 and E 7 ). The 
algorithms for the exceptional groups are difficult to implement and for E& no algorithms are 
given. But if one constructs a program for SU(N) it can be rather fast, and for 'simple' QCD 
calculations this is all one needs, because the quarks are in the fundamental representation 
and the gluons are in the adjoint representation. Such programs can be quite short in any 
symbolic manipulation language (see for instance ref. [|[). For most calculations thus far 
this is all one needs because for QED and SU(2)xU(l) the group theory is completely trivial. 
Nowadays however there is much activity concerning what lies beyond the standard model, 
and a much larger variety of groups can occur. Hence it would be wise to present perturbative 
calculations in such a way that the group and the representation(s) have not been fixed yet. 
Yet one would like to have a compact result. The essence of such an endeavor is or course 
to utilize no group or representation specific information, and to express the result in terms 
of as few invariants as possible. These invariants then can be either tabulated for different 
groups and representations, or there should be easy algorithms to evaluated them. What we 
are going to show here is a very short version of a recent paper Q on this subject. For more 
details the reader should consult the original paper. 
First we take 

[T a R , T b R ] = i f abc T R (1) 
Of special interest are two quadratic Casimir operators: 

i T R T R)ij = CrSij (2) 
jacdjbcd = ^ab (3) 

The index R in Tr and Cr labels the representation. The next set of invariants we consider 
are symmetrized traces. 

Str T ai . . . T an = — Tr T a *W . . . T a ^ (4) 

For each representation one may define a symmetric invariant tensor cIr with 

= Str T* 1 ...Tg (5) 

This does overparametrize the problem, but if we manage to express all color factors in terms 
of contractions of such objects we have made an enormous simplification. 

For a reduction into invariants we first deal with all generators that are not in the ad- 
joint representation. If there are still open indices we multiply with appropriate projection 
operators. This gives complete traces and we have to find an expression for them in terms of 
symmetrized traces. Of course one can eliminate contracted indices inside the same trace as 
one does with traces over gamma matrices. Here this is more complicated though. Writing 
the expression in terms of symmetrized traces can either be done by recursion or with a closed 
formula which is a bit messy to write here. The formula gives much faster results, but in both 
cases the expressions become rather long when there are many generators in the trace. 
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Next one has an expression with invariants d,R and 'structure constants' / which can 
be considered as proportional to generators of the adjoint representation. These do not 
necessarily occur in loops. If they do, they can be written also in terms of invariants d A - 
After this step we are left with combinations of invariants and structure constants in which 
the structure constants cannot be arranged in terms of loops. At this point we start using 
Jacobi identities. For a computer this is not easy, but the program we constructed can do this 
for color traces that contain (at the beginning) up to 14 generators without any problem and 
for 16 it works almost always. This manages to reduce all such traces to contractions between 
invariants d with the exception of one combination of three tensors cIr and two structure 
constants / when we started with 14 generators, and three similar objects at the 16 generator 
level. In some cases these objects can be reduced (like if at least two of the three invariants 
d belong to the adjoint representation, or for groups for which some identities hold which 
includes almost, but unfortunately not all groups). The resulting objects are considered as 
fundamental by the program. It is possible to express them in terms of an even smaller set 
of independent objects, but unfortunately this will not make the expressions shorter, because 
the constants in such a reduction are not simple. Hence the program gives its answers in these 
objects and one can evaluate these contractions afterwards for any given group. The paper 
presents also the formalism on how to do this for each group and for any given representation. 
Here we just give some examples of some rather complicated color traces. All examples use 
an experimental version of the program FORM which will be released later this year. 

We first look at the following trace: 

n nn — 1 ' " ' 1 R 1 R " ' 1 R \ 

which gives some type of maximal complexity. For n = 7 we obtain: 
R 77 = 112/3 d a R bcdef d a A bc9 d d A ef9 - 328/9 d abcdef d abc9 d d A ef9 
+d abcdet 'd abcdet (-56 C R + 296/3 C A ) 

+d a R bcd d A ef d A ef (42 C R - 749/10 C A ) + 67/15 I 2 (R)d a A bcd d a A bef d A ef 

+d a R bcd d a A bcd (35 C% - 357/2C 2 R C A + 868/3 C R C\ - 2695/18 C\) 

+h{R)d a A bcd d a A bcd {7 C\ - 1603/60CrC a + 497/20 C\) 

+N A h{R){+Cl - 21/2 C\C A + 175/4 C R C\ - 280/3 C R C\ 

+5215/48 C 2 R C\ - 19075/288 C R C\ + 43357/2592 C\) (6) 

in which Na is the dimension of the adjoint representation and I 2 {R) is the second index of the 
representation R and it can also be written as I 2 {R) = (NrCr)/Na with Nr the dimension of 
the representation R. This computation took less than 35 sec on a PP200 running NeXTStep. 

If the representation R in the above example is taken to be the adjoint representation 
things are much simpler and much quicker: 

A 77 = -^d a A bcdef d A bc9 d d A ef9 + ^C A d a A bcd d a A bef d cdef - ^N A C A (7) 

and the computer time needed is less than 1 sec. Finally a topologically complicated example: 
It is called the Coxeter graph. It contains 14 vertices and the smallest loop in it has 6 vertices. 

G 6 (n = 14) = ^-d abcdef d A bc9 d d A ef9 - ^C A d a A bcd d abef d cdef + ^N A C\ (8) 
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This took 1.6 sec. 



3. Sums 



The method of evaluating structure functions in Mellin space dates back to the origins 
of QCD ||] It has been used || to obtain the anomalous dimensions of the deep inelastic 
structure functions at the two loop level. Kazakov and Kotikov [Q] used it for obtaining the 
ratio R = ctl/&t of deep inelastic structure functions at the two loop level. A complete 
program for all two loop calculations has however not been constructed thus far. Hence there 
is a challenge here. The main prize would of course be a program that can compute the 
three loop anomalous dimensions. We will address here a project in which first a program is 
constructed to evaluate all relevant two loop coefficient functions in Mellin space. Then one 
has to look at the inverse Mellin transform to obtain results in x-space. Next would be the 
study of three loop anomalous dimensions. Thus far we cannot say much of this. 

Let us have a look at a typical two loop diagram. 

= j d D Pl ^P2 (p?)a((p+pi)2)4 ^ )6((p+p2)2)B(pi)c(p 2 )d(p 2 )e ( 9 ) 

We can attack this diagram in a variety of ways. The method one might prefer is the 'brute 
force' method. One decides to compute the iV-th moment in which we keep N symbolic. 
Therefore the two denominators are expanded, which results in a single symbolic sum. The 
resulting integral can be attacked with the standard techniques, but one has to introduce 4 
more sums. The only good news is that no individual term has more than 4 nested sums. 
Amazingly enough these sums can be solved, although one has to do quite some work teaching 
the computer summation. The summation packages of the big computer algebra programs 
are almost useless here. One runs however immediately into trouble when trying to solve a 




similar toplogy 




Thus far the brute force method has failed on this diagram. 



Hence one has to be a bit smarter. If one looks in the Kazakov and Kotikov paper one can see 
a number of reduction schemes by which relations between the various topologies are derived. 
Two of them can then be marked as the simplest ones. These are then evaluated. They need 
only a single sum if enough preperatory work is done. Then a next level of topologies can be 
done which gives sums that involve these simpler topologies. But because the whole scheme 
is properly built up, the sums do not mix in such a way that we run into very complicated 
sums. Still one may need a number of sums that are not readily available and hence quite 
some attention has to go into the construction of a program that can handle all available 
sums. Because the complete program is still under construction I cannot show too many 
details here but yet, some may be interesting. 

The class of functions that we run into is called 'harmonic series'. What is shown here 
about them can be found in a more complete version in ref Q. There exists a variety of 
notations for them. Because the more conventional notation is not very useful for computer 
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programs we use a slight modification of this notation. The basic function is 

n 1 

<Mn)=£— m>0 (10) 



i 1 
i=i 



n (—] V 

= Y±-L m<0 (11) 

Z-^/ Am v ' 

i=l 1 



and higher functions are defined by recursion: 



n 1 

5 mii ..., mfc (n) = S m2: ..., mk {i) mi>0 (12) 



Ei^r 5 ««.- ) m fc (i) ^i<o (13) 



=1 



These functions appear, amoung others, when T-functions are expanded in terms of e. But 
they also pop up in sums of the type 

n , 

i=i 

which are rather common in these calculations. Certain classes of sums can be evaluated to 
any complexity of the participating harmonic series. An example is combinations of S...(n—i), 
S...(i) and powers of 1/i as in 

n—l ^ 

Y] ^S 2 (n-i)S- 2 -i(i) = -352 -4,-i (n) - 45 2 _ 3 _ 2 (n) + 45 2 _ 3 ,-i,i(n) - 35 2 _ 2 _ 3 (n) 
«=l 

+2^2 -2,-2,1 (n) + 5 2 _ 2 _i i2 (n) + 2S 2 ,2,-2,-i(™) + 25 2 ,3,-i,-i(n) 

-45 3 _ 3 _i(ra) - 4S 3 _ 2 ,-2(rc) + 45 3 _ 2 -i,i(n) 

+45 3 ,i,- 2 ,-i(n) + 25 3 , 2 ,-i,-i(n) - 35 4 ,_ 2 ,-i(n) (15) 

and so on. There are also complicated sums that have thus far resisted generalization. This 
means that we cannot put in arbitrary harmonic series. An example is 

^ZPMr Sl{n+J) = (^iK^2) (25l( ^ (ri - 2) - 52(n - 2) 

-5_ 2 (n+2) - 5_ 2 (n-2) + (_l)»( n 2 +n+3) 

-5i(n-2)(l+— ^- + --— ) (16) 
1 A n-1 n n+1 n+2 J ' v ; 

The sums become also much more difficult when only a 'partial range' is considered as in 

m -. 

EH)'(^ = ??? ( 17 ) 

i=l 
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One solution would be the introduction of a new type of functions with two variables. One 
would then construct the proper relations for these functions and hopefully, at the end of the 
complete calculation most or all would cancel. This approach may be necessary for the three 
loop anomalous dimension. At the two loop level it is not needed. 

The sums that occur in the evaluation of diagrams in Mellin space are of course closely 
related to the sums that one runs into when making Mellin transforms of complete results 
in x-space. However, in the case of the transforms there is an extra class of sums: sums to 
infinity. These are in principle easier, but we need to derive some extra relations for them if 
we would like to be able to do them to a sufficient depth. Additionally one needs to know 
the values of the harmonic series at infinity. These values give us a number of constants that 
do not enter the problem if we calculate the diagrams directly in Mellin space. Hence their 
cancellation serves as a good check. For instance, the Mellin transform of C2 in ref [|| results 
in a formula with 154 terms. Because this is a new result it is shown in the appendix. It does 
take the program however only a few seconds to evaluate it. There is however an interesting 
spin-off. From the calculations in Mellin space, we know which classes of functions can occur 
in the Mellin transform. We can construct a basis in the space of these functions. If we 
have an equal number of functions in x-space of which the Mellin transforms span the space 
formed by this basis, we can do the inverse Mellin transform by just solving a linear set of 
equations. Additionally it tells us that there cannot be any relations between the functions 
in x-space that we used. Let us have a look at this basis. First we define the 'level' of a term 
that involves an harmonic series. This level is the sum of the absolute value of its indices 
and to that we add the number of powers of denominators in the term. Hence the level of 
the argument of the sum in eq. 15 is 2 + 2 + 3 = 7. This is also the level of the terms on the 
right hand side of the equation. For two loop calculations we will need only functions of a 
level up to four. At level one there are only two functions: Sx(n) and S-i(n). The function 
1/n can be written as S\ (n) — Si(n — 1). The two terms in this expression correspond to the 
same function in x-space except for an overall factor 1/x for the second term. Hence we will 
not consider 1/n as an independent function. Additionally we do not have to worry about 
products of harmonic series with identical arguments. These can always be expressed in terms 
of sums of terms that have only single harmonic series as is shown in an example: 



At any given level k greater than 1 there are three times as many functions than at the 
previous level. For each function at level k — 1 one can construct a function at level k by 
adding at the left an index 1, an index -1, or by raising the leftmost index by one (in absolute 
value). Hence at level 4 there are 54 independent functions. Because we also need the lower 
functions one has to consider 80 functions in x-space before one can make a full inverse Mellin 
transform. An example is: 



Si t2 (n) S , _i,i(n) 



-S-2,i,2{n) - 5_ 2 ,2,l(n) + S- 2 ,3{n) + 2 S-i,i,i,2(n) 
+£-1,1,2,1 («) - £-1,1,3(71) - £_i,2,2(n) - Si ,-3,1 (n) 
+£i,-i,i,2(«) + Si ,-1,2,1 (n) - Si -i,z{n) + £1,2,-1,1 (n) 



(18) 



S-\,\,v{n-l) 



Li 3 (l-x) ln(l-x) 1 



1+x 1+x 2 1+x 



+8{l-x){-\u + l -C, ln(2) - ^(ln(2)) 4 - Li 4 (i) ) 



(19) 
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There is one case in which the basis of single higher harmonic series is not practical. This 
is when one has to evaluate these series at infinity. In that case some of these objects can 
be infinite and one would like to cancel the infinities between the various terms. In principle 
all divergent objects can be expressed in terms of powers of just a single divergence Si(oo) 
times finite terms. This is a rather soft divergence which can be regularized rather easily by 
replacing the infinity temporarily by a large integer N. In some case one has to worry then 
about whether objects go to infinity like N or like 2N in which case one gets additional finite 
contributions as in Si(2N) — > S\(N) + In 2. These things are rather straightforward though. 

A. Two loop moments 

The Mellin transform of the coefficient functions c 2 from Zijlstra and van Neerven || 

c 2 = +6(N-3) Si _ 2 (iV-3) (8/5C F C A - 16/5Cf ) 
+6(N-3) 5i_ 2 (iV-2) {-8/hC F C A + 16/5Cf) 
+5(N-2) ( 3 (12/5C F C A - 24/5C 2 F ) 
+9{N-2) (+5i(iV-2) (8/bC F C A - 16/5C|) 
+5 2 (iV-2) (8/5C F C A - 16/5Cf ) 

+S^{N-1) {\2C F C A - 24C 2 F ) + S_ 3 ,i(iV-l) {-8C F C A + 16Cf ) 
+S- 2 (N-1) (8C F C A - 16C|) + S_ 2 ,_ 2 (iV-l) (-24C F C A +48C 2 F ) 
+Si(N-l) (I585/54C F C A - 89/27C F n f + 5/2C|) 
+5i(iV-I) Ca (S6C F C A + 48C F ) 
+Si_ 3 (iV-l) (~24C F C A + 48C 2 F ) 

+5i,_ 2 (iV-l) (3GC F C A - 72C 2 F ) + Si,_ 2i i(JV-l) (8C F C A - 16Cf ) 

+iSi,i(JV-l) {31l/9C F C A - 26/9C F n f - 43C|) 

+5i,i,i(iV-I) (22/3C F C A - 4/3C F n f + 8C 2 F ) 

+S 1 , 1 ,_ 2 (iV_i) (24C F C A - 48C 2 F ) + S hlxl (N-l) (24C 2 F ) 

+S h2 (N-l) {-22/3C F C A + 4/3C F n f - 4C 2 F ) 

+5i,i, 2 (iV-l) (4C F C A - 32C 2 F ) + 5i, 2 ,i(iV-I) (-4C F C A - 24C 2 F ) 

+S 2 {N-l) (-212/5C F C A + 4C F n f + 189/5Cf) 

+5i, 3 (iV-l) (12C F C A + 4C F ) + ,S 2 ,_ 2 (iV-l) (-8C F C A + 16C 2 F ) 

+S 2 ,i(N-l) {-44/3C F C A + 8/3C F n f + 8C 2 F ) 

+S 2AA {N-l) (-24C 2 F ) + S 2 , 2 (N-1) (20C|) 

+S 3 (iV-I) (55/3C F C A - 10/3(7^/ - 18C|) 

+S 3)1 (N-1) (8C F C A + 8C F ) + 5 4 (iV-l) (-12C F C A + 14C|) 

+Si(N) (-4639/45C F C A + 110/9C7pn/ + 337/5C F ) 

+5i(iV) Cs (72C F C A - 144C F ) + S h . 3 (N) (24C F C A - 48C 2 F ) 

+5i - 2 (iV) (-56C F C A + 112Cf.) + _ 2 (iV) (-48C F C A + 96C F ) 



+Si,i(iV) (-68C F C A + 4C F n f + 84C F )+S ls ,i(N) {-8C 2 F ) 
+S lt3 (N) {-2AC F C A + A8C 2 F ) + S 2 ,_ 2 (iV) {-IQC F C A + 32C|) 
+S 2 {N) (7AC F C A - 4C F n f - 7AC 2 F ) + S 2>1 (N) (16C|) 
+5 3 (iV) (-20C F C A + 28C|) + S 3 ,i(N) {-8C F C A + 16C|) 
+S 4 (AO (12C F C A ~ 2AC 2 F ) 

+5i(JV+l) (3914/27C F C A - 488/27C F n f - 121Cf) 
+Si(N+l) Cs (-84C F CU + 144C|) 
+5i,_ 3 (iV+l) (-40C F C A + 80C£) 

+5i _ 2 (AT+l) (20C F CU - 40Cf ) + 5i _ 2j i(iV+l) (8C F C A - 16C F ) 

+5i,i(JV+l) (668/9C F C A - 68/9C F n/ - 68C F ) 

+51,1,1 (JV+1) (22/3C F C A - 4/3C F n/ + 36C F ) 

+Si,i,- 2 (iV+l) (56C F C A - 1126$) + (24C F ) 

+5i, 2 (iV+l) (-22/3C F C A + 4/3C F n f - 32C 2 F ) 

+Si,i )2 (N+l) {AC F C A - 32C 2 F ) + 5i, 2 ,i(iV+l) {-AC F C A - 24C 2 F ) 

+5i, 3 (iV+l) (28C F C A - 28C$) 

+S 2 {N+1) (-±909/15C F CU + 38/3C F n f + 646/5C F ) 

+S 2 ,i(N+l) (-U/3C F C A + 8/3C F n f - A8C 2 F ) 

+5 2 ,i,i(iV+l) {-32C 2 F ) + S 2>2 {N+l) {28C 2 F ) 

+S 3 (N+1) (115/3C F C A - 10/3CW/ - AC F ) 

+53,1 (JV+1) (8C F C A + 2AC 2 F ) + S 4 (^+l) (-12C F C A - 6C 2 F ) 

+(-Si(N+2) + 5i,_ 2 (iV+2)) (72/5C F C A - UA/5C 2 F ) 

+(S 2 (N+2) + S 3 (N+2)) (72 /5C F C A - 144/5Cf ) 

+5i,_ 2 (iV+3) (-72/5C F C A + 144/5Cf ) 

+S 3 (N+3) {-72 /hC F C A + 144/5C F ) 

-5465/72C F C A + 457/36C F n/ + 331/8C£ 

+C 3 (54C F CU - 72C F )) 
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